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Abstract 

Topologically-ordered matter is a novel quantum state of matter observed only in a 



B 

P3 ■ small number of physical systems, notably two-dimensional electron systems exhibiting 

O . fractional quantum Hall effects (7). It was recently proposed (2-4) that a simple form of 

topological matter may be created in interacting systems of rotating ultra-cold atoms. We 
describe ensemble measurements on small, rotating clusters of interacting bosonic atoms, 
j> ■ demonstrating that they can be induced into quantum ground states closely analogous to 

r^ ! topological states of electronic systems. We report measurements of inter-particle correla- 

tions and momentum distributions of Bose gases in the fractional quantum Hall limit, mak- 
^ ■ ing comparison to a full numerical simulation. The novel experimental apparatus necessary 

j-^ . to produce and measure properties of these deeply entangled quantum states is described. 

o 
o 

The current understanding of many-body phases and phase transitions is dominated by the 
concept of spontaneously broken symmetry, as described by a locally-defined order parameter. 



X 

c3 ■ Topological matter is a fundamentally different type of ordered matter that cannot be charac- 

terized by a local order parameter and its long-range correlations (5). Nor can transitions into 
these phases be described by Landau symmetry breaking, which until recently, was thought to 
describe all phase transitions. In the quantum (zero-temperature) phase transition into these 
novel states of matter, particles form deeply entangled, long range states which can have un- 



usual properties, such as excitations with fractional charge and fractional quantum statistics. 
Some topological phases may also have applications in topological quantum computation (6). 

In the fractional quantum Hall (FQH) effect, interacting electrons transition from a nearly 
free electron gas into a new strongly correlated state describable by composite particles of elec- 
trons associated with a number of magnetic flux quanta. An analogous situation has been pre- 
dicted to occur for an interacting Bose gas rotating a frequency Q that approaches the harmonic 
oscillator frequency u of the atoms confined in a trap. In this case, the quantized single-particle 
energy levels become highly degenerate, and the Hamiltonian of the system has a form similar 
to the interacting electron gas that produces fractional quantum Hall states. The lowest en- 
ergy states transition from a weakly interacting Bose gas to highly entangled topological states 
at nonzero and quantized angular momentum. Additionally, the excitations of this system are 
predicted to possess fractional statistical character (7). 

The atomic equivalent of a FQH state may be achieved only when the number of particles 
A^ becomes comparable to the number of vortices A^^ in the ground state, i.e., when the filling 
factor N/Ny is of order one. Long-lived vortices have been excited in superfluid atomic gases 
in a variety of ways (8-12). In previous experiments probing high angular momentum Bose 
condensates consisting on the order of 10^ particles or more (13-15), the lowest Landau level 
was achieved, and a softening of the Abrikosov lattice of superfluid vortices was observed 
(75). However, in these experiments the FQH regime was not achieved, since the filling factor 
was of order 100, and the energetic gap to excited states, scaling inversely with A^ (16), was 
significantly lower than achievable temperatures. 

In this article, we describe the realization of the FQH regime by probing an ensemble of 
interacting, rotating Bose gases in the few-body limit of 1 < A^ < 10. The key development 
in this work is the production of an optical lattice potential, each of whose sites exhibits a 
locally rotating deformation to its nominally isotropic and harmonic confining potential (see 




Figure 1 : Simulation of a magnetic field at the rotating minima of an optical lattice potential. 
(A) Bose-condensed atoms are loaded into a deep optical lattice formed by the superposition 
of a two-dimensional lattice (blue beams) and a standing wave along its normal (red beams). 
The optical phases of the 2D lattice beams are modulated in time to produce lattice sites which 
near the lattice potential minima rotate at a rate comparable to the trapping frequency for atoms 
confined there. (B) The potential near each site approximates a strongly oblate harmonic oscil- 
lator, whose principle axes rotate in the plane normal to the direction of tightest confinement 
at a frequency Q ^ caj± (C) Though the local potential near each site rotates, the overall lattice 
structure remains constant (77). 



Fig. 1). An adiabatic sequence of trap deformation strength and rotation rate is used to attempt 
to populate FQH ground states near the centrifugal limit. Following this manipulation, pulses 
of laser light which photo-associate pairs of atoms into electronically excited molecular states 
are used as a probe of inter-atomic correlation. Additionally, time-of-flight images are taken 
to measure the momentum distribution of particles following the rotation sequence. Both are 
shown to be in qualitative agreement with direct simulation of the few-body system, and provide 
evidence for strong correlations. 

For an interacting gas of bosonic atoms confined to an oblate cylindrically- symmetric har- 
monic oscillator (see Fig. IB, for uj± = uj), the total energy in the rotating frame and the lowest 
Landau level may be written in units of the vibrational quanta Tiuj as {18) 



{mi} 

where Lq = ^^ mal^am is the total angular momentum along the axis of rotation in 
units of h, and a]^ creates a particle in an angular momentum eigenstate m with wavefunc- 



tion (/)m[^= {x + iy)/ao] = z^exp(— |zp/2). Here, ao = Jh/rriaUJ, with rria the atomic mass. 
Excitation along the axis of rotation is effectively frozen out by the condition that the trapping 
frequency in that direction is much larger than the thermal and interaction energies. Interactions 
are assumed repulsive and characterized by the dimensionless parameter r] = a^/v^an, where 



as is the s-wave scattering length for the constituent atoms, a^ = Jh/rriaUJz is the confinement 
length along the rotation axis. The quantity v^^i} characterizes the rate at which two atoms in 
angular momentum states mi^2 are scattered into 7713^4 by interparticle interactions. Ground state 
correlation is expected to enter in the regime ^/2^i{N — l)r] ^ (1 — Q/uo), where interactions are 
sufficient to begin mixing the single particle degeneracy near the centrifugal limit, and become 
strong in the limit ^/27^r] ^ (1 — Q/uj). An energy spectrum for four particles is shown in 
Figure 2; the reduced interaction energy of strongly correlated states in general introduces level 
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crossings between the zero angular momentum ground state and correlated states of higher L^ 
at rotation rates Q < uo. 

In this experiment, small numbers of atoms are first isolated and evaporatively cooled into 
the non-rotating ground states of many optical microtraps formed near the antinodes of an 
optical lattice. Angular momentum is then transferred into the gases by applying a rotating 
deformation to the lattice potential near its antinodes, such that the local potential approx- 
imates a rotating anisotropic harmonic oscillator. This perturbation, describable as adding 
He = Z)m ^ {^m+2^m + Sm^m+2) to the cucrgy Hrot abovc, mixes states which differ by two 
units of angular momentum, and produces avoided crossings in the energy spectrum for the 
few-body states (Fig. 2b). By controlling both the rotation rate Q and deformation amplitude 
6 = 2(cl;+ — cl;_)/(cl;+ + cl;_), an adiabatic pathway is employed to transfer atom clusters from 
zero angular momentum into correlated states of higher angular momenta. It is important to 
stress that rotation is applied only locally to each lattice site, and is not applied to the global 
lattice structure as presented elsewhere (79), allowing application of the rotating trap idea to 
small clusters of atoms. 

At each ground state crossing, in the ideal case of a perfectly harmonic microtrap, Nh 
units of angular momentum enter into a sample consisting of A^ atoms, and simultaneously 
cause the second-order correlation of the gas to decrease by a significant degree. This can be 
understood as the few-body equivalent of successive vortex nucleations, each accompanied by a 
strong redistribution of the few-body wavefunction away from the mean-field to accommodate 
energetic constraints of repulsion. This continues as a sequence of tabulated (3, 7, 16) ground- 
state crossings terminating at a total angular momentum of A^(A^ — l)/i, where the final state 
closely resembles (3) the |-Laughlin (j){zi) = Ylij{zi — ZjY exp(— Y.^ l^iP), a bosonic analog 
of states (20) previously discussed in context of electronic systems, and where the probability 
to find two particles in the same location is zero. 
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Figure 2: Energy spectrum and excitation gap as a function of trap rotation rate Q and de- 
formation strength e for four particles. (A) Energy spectrum for four particles in the lowest 
Landau level of a rotating harmonic trap. Without interactions (left, r] = 0), states collapse into 
a macroscopic degeneracy at the centrifugal limit Q = uo. States with total angular momen- 
tum in integer multiples of the particle number are shown as solid lines L = (red),4 (blue), 8 
(green), 12 (black). Repulsive interactions (middle, characterized by dimensionless parameter 
T] = 0.0054) shift the energy of uncorrelated states upward, perturb degeneracy, and intro- 
duce level crossings in the ground state contour toward increasingly correlated states. Adding 
a rotating quadrupole deformation (right, e = 0.016) to the harmonic potential couples states 
whose total angular momentum differ by 2h and produces avoided level crossings. The ground 
state contour provides an adiabatic pathway into correlated states near the centrifugal limit with 
nonzero total angular momentum. (B) shows a small scale view of the the first ground-state 
crossing for four particles from a trivial non-rotating state to the four-particle Pfaffian with 
L = 4. (C) shows the excitation gap as a function of both Q and e. The full adiabatic rotation 
sequence is illustrated as a dotted line, terminating at a final frequency f^/; the entire sequence 
is translated in frequency as a control parameter. 



As a specific example, for A^ = 4 particles as shown in Fig. 2B, the ground- state transitions 
first from an uncorrelated L = state (p{zi) = Hi exp(— |z^p), through a state of L = 4 previ- 
ously identified (3, 16) with a Pfaffian-type wavefunction (27) at f] = 0.932cl;. The subsequent 
transition is to an L = 8, |-quasiparticle (76) state at f] = 0.9776l; before the final L = 12, 
i-Laughlin state at f] = 0.9896J. 

Atoms that adiabatically cross into the fractional quantum Hall ground states are expected to 
show a decrease in second-order correlation. The degree of correlation in the few-body atomic 
states is probed with photoassociation to electronically excited molecules (22) (see methods). 
The molecular formation rate F oc / n^{z)g2{z)d\, where the second-order correlation function 
g2[z) = lim^/^^(?/i|V^^,4^4)/(^|4)(^l/4^) (23), with 4 = Em0m(^)Sm the atomic field 
operator, measures the likelihood of pairs of atoms to be found at short range, compared to that 
expected from a mean-field form at the same density. The number of atoms surviving such a 
photoassociation pulse, therefore, directly probes the degree of correlation in the final few-body 
state of the atoms. 

The optical lattice employed here (see Fig 1 A and Methods) is the superposition of a two- 
dimensional triangular lattice potential formed by the intersection of three laser beams of equal 
frequency and a separate one-dimensional lattice formed by a separate counter-propagating 
beam pair along the rotation axis. The local lattice site structure approximates a strongly oblate, 
cylindrically symmetric three-dimensional harmonic oscillator potential with trapping frequen- 
cies of {u)±,ujz) = 27r X (2.1, 28)kHz. This results in relatively strong interactions between 
atoms, with interaction parameter rj = 0.0054, separately verified by measurement of collec- 
tive mode frequencies (see Methods). Care was taken to produce nominally symmetric lattice 
sites with uo^ c^ uj^ to an accuracy of 0.5% over a day. The method used to generate a con- 
trolled rotating trap deformation consists of diabatic translation of the two-dimensional lattice 
potential obtained by rapid electrooptic modulation of lattice beam phases, and is described in 



detail in the methods section below. The time-average lattice potential has a static overall lat- 
tice structure, whose sites approximate locally harmonic oscillators with dynamically controlled 
anisotropy e and orientation of their principle axes. 

Atoms are loaded into this potential from a magnetically trapped Bose-condensate of ^^Rb 
atoms, by first adiabatically increasing the two-dimensional lattice potential to confine atoms 
into tubes formed by the intersecting blue beams shown in Fig. lA, then allowing the gas 
to expand along the axial direction. The axial standing wave formed by the red beams (Fig. 
lA) is then adiabatically increased to a final depth sufficient to prevent tunnelling along the 
axial direction for relevant timescales. This is followed by state- selective evaporation of atoms 
in lattice sites at the outer edges of the initial Bose Einstein condensate (see Methods), and 
selective evaporative removal of atoms to a final mean density of A^ ^ 5 atoms/site in the full 
three dimensional potential. This results in a total number of atoms of order 10^, distributed 
over several thousand lattice sites. 

The local lattice potential is then rotated with a time-dependent amplitude and frequency 
ramp determined to be a sufficiently adiabatic pathway to correlated states using a full numerical 
diagonalization of the hamiltonian in eq. 1, as described in the Methods section. The few-body 
energy spectrum and gap to first excited state are presented in figure [21 for A^ = 4 atoms, along 
with the chosen ramp sequence (Fig. 2C). The total ramp time used in this experiment is 48ms, 
consisting of three 16ms piecewise linear portions connecting the points {{Q. — Q.f)/uj,e) = 
(-0.065, 0), (-0.025, 0.06), (+0.005, 0.04), (0, 0). The sequence was translated in frequency, 
to allow the final rotation rate Vtf to be used as a control parameter. We note that the adiabatic 
ramp should serve to couple multiple occupancy classes into correlated states simultaneously, 
as each level of occupancy A^ carries a similar, but different ground state crossing sequence. 

The final state is probed by simply releasing the atoms from all applied potentials. The 
resulting atomic distribution is then imaged along the rotation axis after a 10ms time-of-flight. 
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Figure 3: Angular momentum distribution of atoms in rotating lattice potential. Absorbtion 
images (A) probe the angular momentum distribution following adiabatic rotation sequences at 
variable final rotation rates. As the rotation rate is increased, the cloud profile first grows slightly 
in radius, followed by the development of a halo of high angular momentum components. Close 
to the centrifugal limit, the distribution broadens rapidly. Plotting density as a function of radius 
from the center-of-mass (B), the momentum distribution can be fit well to the form described 
in the text. (C) The expansion coefficients c^ display a nontrivial behavior as the centrifugal 
limit is approached from below. The dependence of the lowest angular momentum components 
Co and ci + C2 are compared to numeric predictions (blue lines). In the blue region, marking 
the range spanned by the first ground-state crossing location for particle numbers 2 < A^ < 6 
(vertical lines), we find reasonable quantitative agreement with the model, implying successful 
generation of FQH states of angular momentum A^. The small discrepancy seen at rotation 
rates Qf ^ 0.88 are expected from a truncation of numeric predictions at occupancies N > 7; 
incorporating a phenomenological prediction (dotted curve, see Methods) for A^ = 7 improves 
the fit. 



Due to the tight confinement of the lattice, the evolution in time-of-flight is primarily determined 
by the initial momentum distribution before release and is largely unperturbed by interactions. 
Furthermore, due to extremely weak tunnel coupling between sites, the time-of-flight density is 
an incoherent sum of the momentum distributions in all sites. In the lowest Landau level, time- 
of-flight expansion proceeds such that the momentum distribution mimics the in-situ density 
n{z) with \z\ replaced by p/po (24), where p is the magnitude of momentum in the plane of 
rotation, and po = {mahuoY^'^ is the oscillator momentum. This leads to a momentum space 
distribution of the form n{p) = ^^ Cmip/poY^ exp{—p'^ /pl)/m\, where c^ for m = 0, 1, 2... 
characterize the angular momentum distribution. 

At sufficiently low final rotation rates, the adiabatic pathway terminates prior to the first 
ground state crossing, and the momentum distribution should reflect the zero angular momen- 
tum ground state. For rotation rates Qf < 0.87cl;, we observe a gaussian time-of-flight profile 
of momentum width rria x 3.8mm/s, consistent with atoms remaining in the ground state, given 
Po ^ ^a X 3.1mm/s, the imaging resolution and initial spatial extent of atoms across the lat- 
tice. This verifies an ability to provide controlled rotation sequences without unanticipated 
heating. The angular momentum distribution c^ can be extracted from images using singular 
value decomposition to the fit form for n{p) above, and as described in the methods; the result 
is shown in figure 3C. At rotation rates 0.87 < ^f/oo < 0.94, the gaussian distribution broad- 
ens, indicating increased occupancy of single-particle states of angular momentum m = 1. For 
Qf > 0.95cl;, a high momentum "halo" develops with significant population in m > 6, growing 
in amplitude for frequencies 0.950 < ^f/oo < 0.980. These halos are unexpected in the ide- 
alized case of a perfectly harmonic oscillator potential, but appear in the coherent evolution of 
higher occupancy clouds for an anharmonic rotating trap with a negative quartic component, as 
described in the methods below. Beyond a final rotation rate Qf > LOlcj, the central peak in 
time-of-flight grows in size rapidly, reflected by a growth in components c^ for 1 > m < 4. 
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The expected form of nN{p) was calculated from full dynamic evolution of the hamiltonian 
presented above for all occupancies A^ < 7 as described in the Methods. A weighted sum of 
the nN{p) was taken to represent the expected distribution of occupancies in the lattice, and 
is shown in figure 3C, compared to experimentally recorded profiles. In general, one finds 
reasonable quantitative agreement between model and experiment over the range 0.91 < f^/ < 
0.97; the behavior in this range is expected from numerics (see methods) to be dominated by 
the first ground-state level crossing of successively lower A^ systems as Qf increases. Each 
of these crossings is expected to introduce significant correlation, and its precise location is 
determined by the degree to which the repulsive interaction energy is reduced by correlations. 
For A^ = 4 particles, this level crossing has been associated with the entrance of a few-body 
Pfaffian ground state (3, 16). The disagreement at 0.86 < Vtf < 0.90 and slight mismatch in 
slopes for 0.91 < Qf < 0.97 is understandable, due to a truncation of the numeric predictions 
for A^ > 6 where calculations become too numerically intensive; this can be improved assuming 
a phenomenological form for N=7 as described in Figure [3] and methods. The upturn in the data 
very close to the centrifugal limit Qf > 0.97 is currently not understood, but occurs where the 
energy spectrum is exceedingly complicated and sensitively dependent on trap shape. 

To probe atom-atom correlation in the final state, pulses of light tuned to a photoassociation 
transition were applied immediately following the adiabatic rotation sequence. This transfers 
pairs of atoms in close proximity into electronically excited molecules, which quickly relax to 
deeply bound states with sufficient excess energy to be ejected from the trap. The remaining 
atoms were released within 100/iS and imaged in time-of-flight. The total surviving atom num- 
ber Npa then reflects the probability to find atoms in the ground state at short interatomic range. 
Details of this technique appear in the methods and reference (22). 

From Fig. 4, one can see a strong suppression of photoassociation loss near the centrifugal 
limit, consistent with the introduction of correlation in the few-body wavefunction. Here we plot 
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Figure 4: Probing atomic correlation with photoassociation. (A) Shows fraction {N^pa — 
Npa)/Nnpa of atoms lost after rotation sequence and short pulse of photoassociation light. The 
strong suppression at the centrifugal limit is indicative of strong correlation. (B,C) Measure- 
ment of density profile with and without photoassociation allows extraction of the second-order 
correlation function g2{z) as defined in the text. The depression near the centrifugal limit is 
most apparent at nonzero radius. 

the fraction of atoms lost to photoassociation as {N^pa — Npa)/Nnpa, where Npa(Nnpa) is the 
number measured with (without) the photo-association laser pulse. To establish that the change 
in loss rate is due to change in correlation, and not simply decrease in overall density, we form 
a simplistic estimate of the second-order correlation function g2{z, Qf) as defined above. Since 
one expects a constant g2{z, Qf) = {N — l)/N for any A^-particle mean-field state of the form 
(Zlm ^mSm)^|0)' independent of the particular state defined by the 6^, a reduction of g2{z, Qf) 
near the centrifugal limit would provide evidence for departure from the mean-field into the 
fractional quantum Hall regime. We model the loss rate as no{z, t) = —Kg2{z, Qf)nl{z, t), 
where no is the atomic density on a given site, K is a rate constant, and t is the pulse time. 
Using this, an estimate of g2 can be formed as g2{z, Qf) ^ [nQ^{z, t) — n^^(^, 0)]{Kt)~^. The 
quantity Kr is eliminated by forming g2{z) = g2{z, Qf)/g2{0, 0), where 5^2(0, 0) is estimated 
by averaging 5^2(0, Qf) for all points with Qf < 0.875cl;; we find the results plotted in Figs. 4B 
and 4C. 
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In this analysis, we have treated the density after time-of-flight as if it arose from a well 
defined occupancy in the lattice; in reality the time-of-flight density should be treated as an 
incoherent sum over the density nN{p) arising at different occupancies A^, and the estimate 
g2{z) as a type of occupation-averaged correlation. The presence of systems with lower par- 
ticle number contributes a background of uncorrelated low angular momentum states near the 
centrifugal limit, making the reduction in g2{z) there most apparent at nonzero radius. A more 
sophisticated model could give a more realistic interpretation of the quantity 5^2 (^) in terms of a 
weighted average of the second-order correlation of each occupancy class. By systematic varia- 
tion of the average occupancy in the lattice, one could hope to extract separately the correlation 
function for each occupancy. 

In summary, an apparatus has been described to bring a quantum degenerate system of 
bosonic atoms into the quantum-Hall regime by applying a precisely controlled trap rotation. 
We have designed and implemented an adiabatic sequence of trap deformations to populate 
highly correlated states of strongly interacting atoms at nonzero angular momentum, and probed 
its efficiency by measurement of both momentum distributions and interatomic correlation using 
photoassociation techniques. Quantitative agreement between experimental observations and 
full numeric models indicates the establishment of a viable platform to further study entangled 
atomic states in the fractional quantum Hall regime. 
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1 Methods 

1.1 Rotating Lattice Production 

To produce an optical lattice of locally rotating potentials, three laser beams of equal intensity 
and detuned far from atomic resonance are combined with their propagation directions evenly 
distributed on the surface of a cone with a small apex angle. The optical interference pat- 
tern created by these beams consists of a triangular lattice of intensity maxima, whose light 
shift form a conservative trapping potential for atomic motion. Near to the minima, the po- 
tential is locally harmonic, approximately cylindrically symmetric, and may be described by 
l/(x, y) = — Vo Zlj cos(v^/Cerj + 0^), where Tj = cos(27rj''/3)x + 8m{27rj/3)y, with x, y co- 
ordinates in the lattice plane, k^ = 27rsin6^/A, and (pj represent the relative optical phases of 
the three beams (Fig. lA). By choosing a small intersection angle 9 = 8 degrees, the spacing 
between lattice sites is 3.5/xm. This reduces the tunneling rate of atoms between lattice sites in 
the 2D potential to be negligible on experiment timescales, and simultaneously makes the po- 
tential effectively more harmonic by reducing its vibration frequency at a fixed total depth. The 
three lattice beams originate from a common 1.5W, fiber-coupled beam, intensity stablized by 
an acousto-optic modulator and derived from a lOW single-mode Nd:YAG ring laser injection- 
locked to a stable 0.5W non-planar ring oscillator (Lightwave NPRO). To produce a nominally 
cylindrically- symmetric potential near the bottom of each lattice site, center-of-mass vibration 
frequencies are measured with an accuracy of 0.3% along two directions, and beam intensities 
adjusted to equalize these frequencies to a typical precision of 0.5% over a typical experiment 
run time of several hours. The inhomogeneous spread of vibration frequencies due to the gaus- 
sian intensity profile of the lattice beams is calculated to be 0.3% over the loaded volume. 

The locally rotating potential wells are created by inserting two electrooptic phase modula- 
tors into two of the beams forming the 2D lattice potential. Manipulating the relative phases (j)j 



16 



of the beams causes a 2D translation of the lattice potential. By rapidly scanning the potential 
along a given direction (j)j = A(pco8{9 — 27r j/ 3) 8m{ujRFt) at a frequency ujrf ^ ^,^z, 
the lattice potential may be effectively averaged over the timescale for atomic motion cre- 
ating a local potential with a reduced vibration frequency along the axis of translation (de- 
fined by 9). By slowly pivoting this axis in time according io 9 = f^t at a rate of order 
Q ^ 0{u)r), the time-averaged local potential approximates an anisotropic harmonic oscilla- 
tor whose principle axes rotate at rates which can be comparable to the trapping frequency. 
For the data shown here, ujrf = 27r x 500kHz, which was separately verified to be suffi- 
ciently diabatic that heating did not occur on relevant timescales by measuring cloud width 
in time-of-flight. Averaging over the short timescale, the lattice potential can then be written 
as V{x, y) = — Y^j Vj cos(v^/Cerj), where Vj = VoJo{A(j)/2 8'm{Qt + 27rj/3)), and Jq is the 
zeroth-order Bessel function. The local potential near each site minimum is then approximately 
(dropping a constant) V{r, (j)) ~ (mcL;^r^/2)(l + 2ecos(2(/))) in cylindrical coordinates (r, 0), 
with uj = cl;o(1 — A(/)^/32) and e = A(/)^/64, with ujq the radial vibration frequency of a lattice 
site without modulation. The adiabatic pathway described in the text is achieved by controlling 
the amplitude Ac/) of RF phase modulation and the pivot frequency Q.. 

An additional ID potential is used to isolate small clusters of atoms and enhance the effect 
of interactions between atoms loaded into this potential using two additional beams, frequency- 
offset from those forming the 2D potential and counterpropagating along the rotation axis. The 
final potential is a three-dimensional array of highly oblate three dimensional harmonic os- 
cillators, with radial trapping frequencies of up to 6kHz and axial frequencies up to 30kHz. 
Interaction strength is characterized by the ratio of scattering length to oscillator length from 
confinement along the rotation direction, tj = asl\f2j\a^ ^ 0.0054. The parameter t\ was 
separately checked by measuring the variation of the collective quadrupole mode frequency 
with mean particle number {Ml), over the range 3-100 atoms/site, showing a characteristic 
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curve from the hydrodynamic to single-particle regime consistent with the expected interaction 
strength and separately measured center-of-mass frequencies. 

1.2 Preparation of Initial State 

Atoms are loaded from a ^^Rb Bose-Einstein condensate of 10^ particles in the |F = 2, rriF = 2) 
state at a temperature of roughly ZQnK formed by evaporative cooling in a time-orbiting- 
potential (TOP) magnetic trap. After evaporation, the two-dimensional lattice potential is adi- 
abatically increased from zero intensity to its full value of 500mVI^ per beam, after which the 
TOP trap is deformed into a quadrupole trap whose center is pulled below the position of the 
atoms loaded into the tubelike 2D lattice potential. The axial confinement of atoms trapped 
in the two-dimensional lattice potential is adiabatically decreased by reducing the magnetic 
quadrupole field, during which time the axial size of the cloud increases from an a Thomas- 
Fermi radius of 20/xm to a half- width of approximately 200/xm. Following this, the axial stand- 
ing wave intensity is increased slowly to a depth sufficient to inhibit axial motion. In order to 
produce a more defined mean occupancy in the full three-dimensional lattice potential, a tomo- 
graphic technique is used to remove atoms in the tail of the density distribution along the axial 
direction. A weak magnetic field gradient is applied, and a microwave field is applied to transfer 
atoms from the |F = 2, rriF = 2) state into |F = 1, rriF = 1), whose magnetic moment is oppo- 
site in sign. By slowly sweeping the microwave frequency, atoms are adiabatically transferred 
between internal states at the edges of the cloud. Following this, the two-dimensional lattice 
depth is reduced, and a strong magnetic field gradient is applied to completely remove atoms 
in the state |F = 1, ttif = 1) while slowly evaporating atoms in |F = 2, ttif = 2) until the 
desired mean occupancy is reached, as inferred from absorbtion imaging performed transverse 
to the rotation axis. The density profile along the axial direction after this process resembles 
a top-hat, with a physical width of 150/im. Following this, the depth of the two-dimensional 
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lattice is increased to its full value. 

1.3 Analysis of Time-of-Flight 

Time-of-flight momentum distributions were obtained by releasing atoms from all applied po- 
tentials suddenly, allowing free expansion of the atoms for 10ms, then performing absorption 
imaging along the rotation axis. The flight duration was chosen as a balance to obtain a domi- 
nantly momentum- space picture of the cloud, while preserving sufficient density for acceptable 
signal to noise. Profiles were fit to the form n{p) = Z)mCm(l^/Po)^^exp(— p^/Po)/^! for 
the coefficients c^ using a singular value decomposition. The value of po was extracted from 
data taken at small f^/, and fixed for remaining data. A constraint of positivity was enforced 
on the coefficients c^ to remove small negative values by minimizing X]^^/ ScmCmm'Scm^ + 
Y,m -^(cm + ^CjnY with A = 1, whcrc Cmm' IS the covariance matrix obtained from singular 
value decomposition. 

1.4 Photoassociation-based Probe of Short-range Correlation 

Following the adiabatic evolution sequence, correlation is probed by driving pairs of colliding 
atoms into bound, electronically-excited molecules with a short pulse of light tuned from the 
|2, 2) free-atom continuum to the |0^, u = 1, J = 2) molecular state (22, 25) at 12,789 cm"^. 
Light was derived from a external-cavity-diode laser frequency stabilized to ^ 7MHz to an 
optical transition in molecular iodine. After photoassociation, rapid relaxation to the ground 
state is assumed to be accompanied by sufficient kinetic energy that the resultant atoms are 
lost from the lattice with a high probability, and hence the loss rate from the trap during the 
pulse reflects the probability to find two atoms at short range (determined by the outer classical 
turning point of the excited molecular state at tq ^ 2nm (25)). Furthermore, the pulse time 
is made short (100/xs), such that no significant dynamic evolution of the few body state occurs 
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during the exposure to photoassociation light aside from loss. After the photoassociation pulse, 
atoms are released in time-of-flight, and absorptively imaged to obtain both their total number 
and momentum distribution. The time-of-flight image is recorded in successive experiments 
both with and without the photoassociation strobe but the same adiabatic preparation sequence, 
and compared. Total atom number (presented in figure 4A) is extracted by direct integration 
of the raw absorbtion images without additional fitting or modeling. To extract the estimate of 
second-order correlation, the density distribution before and after time-of-flight is decomposed 
according to the fitting method described above. 

1.5 Modeling Adiabatic Evolution 

To model the few-body states, their dynamic evolution, and degree of correlation, direct di- 
agonlization of the few-body hamiltonian Hrot (^q- 1 in the text) including the effects of the 
rotating trap and short-ranged repulsive interactions was performed, similar to that presented 
in Ref. 16. Interactions were modeled by a contact potential, leading to the form of the scat- 
tering matrix v^m^} = C?+m2 (^i + ^2)! 2^1+^2 (fli mi!)"^/^ where 5j is a kroneckar delta. 
Direct diagonalization for fixed e and Q were used to guide the choice of adiabatic evolution 
in trap rotation rate and deformation strength presented in figure 2. The full evolution of the 
few body states, using a time-dependent Shrodinger's equation with varying e(t) and f^(t), was 
used to calculate the expected momentum distribution presented in figure 3C, including effects 
of nonadiabaticity. The effects of trap anharmonicity on the spectrum and dynamic evolution 
were included by modifying the single-particle energy levels to first order according to the ex- 
pected form of the lattice potential at finite depth. 

To account for multiple occupancy classes in the dynamic evolution of the angular momen- 
tum distribution, a separate few-body numeric evolution was performed for each occupancy 
class A^, (results are shown in Figure 5), and a weighted sum of the angular momentum distri- 
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Figure 5: Modeling time-of-flight signatures. (A) The expected distribution P^ of lattice site 
occupancy A^ across the lattice volume with (without) tomographic selection is shown in red 
(black). (B) A separate numeric calculation for the full few-body quantum state is performed 
for each occupancy A^; here the lowest angular momentum component Cm=o is shown for occu- 
pancies 2 < A^ < 6. Each occupancy displays a large drop at the position of the first ground 
state level crossing, followed by a steady redistribution of atoms to higher angular momentum. 
For N > 7, numeric calculation is too intensive, and a phenomenological form is assumed 
(dashed line). Weighting the results in (B) by the distribution in (A), and accounting for effects 
of imaging resolution and cloud size, we obtain the curve presented in Figure 3C of the main 
text. Calculations were performed for the adiabatic ramp parameters described in the text, and 
T] = 0.0054. 
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bution was taken for the results shown in Figure 3C of the main text. To approximate the distri- 
bution of occupancies A^, the number of lattice sites containing A^ atoms was taken to be Poisson 
Pn{^) = ^^ exp{—iy)/N\ with mean number u. While it is in principle possible the preparation 
sequence favors a sub-Poisson distribution through interaction-driven localization, we assume 
this (at finite temperature and over the range of filling factors produced) to be a small effect on 
the full distribution P^. The mean u was assumed to vary spatially across the cloud, following 
a Thomas-Fermi profile in the radial direction, and approximately flat along the axial direction 
due to the tomographic selection described above. For an inverse quadratic density profile of 
maximum filling z/q, the distribution then assumes the form P^ oc Jq^ du u^ exp(— z/)/A^!. The 
peak filling z/q ^ 10 was estimated from absorbtion imaging. The effect of finite imaging res- 
olution and initial cloud size is accounted for by assuming the time-of-flight profile suffers a 
gaussian blur of 1/e-radius IZ ^ 0.7po estimated from cloud profiles measured at low Qf, and 
adjusting the numeric prediction of the density distribution presented in Figure 3C of the text 
according to Cm=o -^ Em CmC^m/ Em ^m, whcrc a-^ = {pq^ + n-^)'^^^pl'^n^. 
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